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Abstract. This paper develops an algorithm to search for approximations to 
a real { which finds rationals of the form a/f)^, such that |^ — < ^^^i in 

B^/^ polynomial time operations, for b < B. An heuristic approach and com- 
putational data give support to the conjecture of O. Karpenkov stating that 
there are infinitely many solutions in integers a and b for such approximation. 



1. Introduction 

Recently, O. Karpenkov :4| posed the problem of the approximation of reals 
by rationals of the form as a natural subject of number theory and for its 
incidence in the solution of cases of the Shrodinger equation. He presents a proved 
lower estimate for the approximation of a real ^ as 



< 



63 In 



and conjectures that there exists a constant c(^) such that the inequality 

c(0 



< 



63 



has infinitely many solutions in integers a and 6. 

In this paper, we study the problem from an heuristic point of view and de- 
velop an algorithm that can compute an approximation of the form a/b^ in Vb 
polynomial time operations. The heuristic approximation to the problem and the 
computational data obtained support the conjecture of O. Karpenkov. 



2. An heuristic approach 

Following a brute force scheme to search for approximations of the form a/b'^ 
to an irrational ^, we can compute, for b — 1, 2, • • • , the optimal value of a as the 
nearest integer to 6^^. For each 6 value we get an approximation 

6_ 
62^ 

It seems reasonable to suppose that the values of 6 are distributed uniformly in this 
interval. The probability of getting a value S < c/b^'^^ , where c is a constant, will 
be, under that supposition, p{b) = 2c/b^^^ . Note that, in this case, 
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Moreover, we can estiniate the number of rationals with the above approximation 

as 

oo oo 
6=1 6=1 

Since the summatory converges for e > 0, we can guess that there may be only 
finitely many rational approximations of the form a/b^ of any real ^, for e > 0. 
When the exponent in h is exactly 3, the summatory is the harmonic series and 

6=1 

where 7 — 0.57721 • • • is the Euler's constant (see [2] Theorem 6.10, p. 136). 
Therefore, the number of approximations with b < B may be estimated as 2c(7 + 
\nB). Note that the fractions of that form are not necessarily irreducible in this 
estimation. 



3. The algorithm 
We consider the following modular quadratic equation, 



(1) 



(mod Q). 



Since Pb^ — aQ + a, for some integer a, we can write 

P a a 

Let ^ be a real number. From the equation 

we have that 
(2) 

If we choose P/Q as a convergent of ^ computed from its continued fraction expan- 
sion, we have that 
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(3) 
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from the well known result due to Hurwitz (see, by example Theorem 193, p. 164 

mm- 

Suppose that we are able to find a solution to equation ((!]) such that b < kQ'^/^ 
and |a| < tQ^I^ for certain constants fc and t. Then we have that 



(4) 



< 
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If we solve equation ^ for some value of a, we have a probability kQ^^/^ of 
getting a value 6 < Q^^^, under the supposition that the quadratic residues are 
distributed at random. Since we check 2Q^/^ values of a, we have the chance of 
getting 2tk values of 6 in average which gives approximations of ^ that fulfill ([4]). For 
this task, we must solve 2tQ^/^ quadratic equations for which there are polynomial 
time algorithms whenever the factorization of Q is known. Suppose that a value of 
6 is obtained such 6 = Q^^^^ so the number of quadratic equations solved to obtain 




it is 0{h^/^). Take into account that there are methods of factorization that can 
factorize Q in Q" steps with u weh below 1/3 (see [6J p. 218), then the complexity 
of this stage is smaller than the solution of the quadratic equations. 

Using the described method, there may be approximations fulfilling ([4]) that 
scape the search of the algorithm. First, a convergent P/Q may have an approxi- 
mation to ^ better than ([3]), mainly for small Q, what allows that higher values of h 
keep bounded the first summand in ([2]). Besides, both summands can compensate 
if they are of opposite sign. For these reasons, in practice, it is better to use a 
brute force scheme for the fist values of h (say for h < 1000) and then to apply the 
described algorithm. Likewise, it is better to check the values of \a\ < Q^^^^^ (by 
example of \a\ < Q"-*^^). Since, once computed the roots of (H]), checking the values 
of b is less time consuming, it is better to rise the bound of 6 to 6 < Q^^^. With 
these trimmings, the algorithm finds all items existing (checked by brute force) at 
least for 6 < 2 * 10^. 

The algorithm was programmed using the PARI package [5] and compiled with 
the utility gp2c. The executable was launched in an AMD64 +3200 computer over 
one week, searching for rational approximations of tt such that the difference is 
less than 1/6"^ in absolute value. There were found 94 rational approximations 
of the prescribed form which can be found in [3] (note that only b values are 
tabulated because the numerator a can be computed from it). The higher rational 
approximation found correspond to 

a = 36266840658555398816245943123914613560. 

b = 3397660065732068041. 

The following figure plots the number of items versus the expected value 1 + 2 In 5 
showing that both agree quite well. Another irrationals as e, \/2, the golden mean 
(1 + •\/5)/2, even numbers suspicious of irrationality as the Euler's constant were 
tested with the same results. 



4. Approximating irrationals by rationale of the form a/p, p prime. 



O. Karpenkov, in his cited paper, also raised the problem of approximating a 
real with a rational with prime denominator. We follow the procedure described 
in the preceding section, to find approximations of that kind. Consider the linear 
modular equation 

(5) Px = a (modQ). 

Since Px — aQ + a 



(6) 
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We take P and Q from a convergent of <^ and solve the equation ([5]) for the lowest 
absolute values of a. Since the density of primes is about 1/lnn, we may expect 
that we will find a x prime for a of the order of biQ and let that prime be p. This 
constructive approach makes reasonable the 

Conjecture. For any real ^, there exists a constant c(^) such that the inequality 

c(0 Inp 



P 



< 



with p prime, has infinitely many solutions. 



Dedicatory: I dedicate this paper to the memory of my friend Manuel Torres 
Hernanz. Since my first collaboration with him, two decades ago, computing the 
Fourier transform of his magnetic measures till a recent talking about the different 
ways to get rational approximations to \/2, I have enjoyed his knowledge, creativity, 
and humanistic spirit. I miss him and I am very indebted to his generosity. 
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